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Abstract 

We calculate, in the framework of open quantum systems, the ground state energy-level shift for 
a static two-level atom outside a spherically symmetric black hole in interaction with fluctuating 
massless scalar fields in the Boulware and Unruh vacuums. We find that the energy-level shift is 
position dependent and thus gives rise to a force on the atom besides the classical gravitational 
force. For the case of the Boulware vacuum that represents a star which has not collapsed through 
its event horizon, this force is attractive near the horizon and is repulsive far away from the black 
hole with a behavior of r~ 3 . For the case of the Unruh vacuum which represents a radiating black 
hole, we find that the contribution to the Casimir-Polder-like force due to the presence of Hawking 
radiation is always attractive and, remarkably, this attractive force diverges at the event horizon. 

PACS numbers: 12.20.Ds, 03.65.Yz, 04.62. +v, 04.70.Dy 
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I. INTRODUCTION 



The shift of the energy-levels of an atom which is caused by the coupling of the atom 
with the quantum vacuum is one of the striking manifestations of the existence of zero-point 
fluctuations, and it is known that this shift can be modified by the presence of cavities jl] 
and the non-inertial motion of the atom itself (2- 5] . One of the observable examples of 
such shifts is the Casimir-Polder force between a neutral electric polarizable atom and a 
conducting plate, which is a result of position-dependent energy-level shift caused by the 
modification of vacuum fluctuations that arise because of reflection of vacuum field modes 
at the conducting boundaries. 

So far, the Casimir-Polder force, which is a consequence of the position-dependent level 
shifts of an atom, is regarded as a result of the reshaping of vacuum fluctuations induced by 
the reflection of field modes at boundaries in flat spacetimes. Since the vacuum field modes 
are also scattered by the curvature in a curved spacetime, such as in the background of a 
black hole, a question naturally arises as to whether the scattering of vacuum field modes off 
the curvature will also produce a Casimir-Polder-like force on an atom outside a back hole, 
and this is exactly what we are going to address in the present paper. 

Our calculation of the energy-level shifts is based upon the framework of open quantum 
systems js], where a static two-level atom outside a Schwarzschild back hole is treated as 
an open quantum system in interaction with the external reservoir of fluctuating vacuum 
massless scalar fields. As for any open system, the full dynamics of the atom can be obtained 
from the complete time evolution describing the total system (atom plus a reservoir of 
external vacuum fields) by integrating over the field degrees of freedom, which are in fact 
not observed. It is worth noting here that the open quantum system approach has also 



been applied to studying other quantum effects in curved spacetime, such as the £ 



awking 
8|- 



radiation of a black hole [7] and the Gibbons-Hawking effect of de Sitter spacetime 

When a curved space-time is considered as opposed to a flat one, a delicate issue arises as 
to how the vacuum state of the quantum fields is specified. In this paper, we deal with two 
vacuum states of the scalar fields; namely, the Boulware vacuum and the Unruh vacuum. Let 
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us note that the Boulware vacuum would be the vacuum state outside a massive spherical 
body of radius only slightly larger than its Schwarzschild radius and the Unruh state that 
best approximates the vacuum following the gravitational collapse of a massive body to a 
black hole. So, we are going to compute the force on a static atom as a result of the modified 
vacuum fluctuations both outside a star which has not collapsed through its event horizon 
and a black hole. 

The paper is organized as follows: In the next section, we will review the basic formalism 
of open quantum systems, the derivation of the master equation describing the system of 
the atom plus external vacuum scalar fields in the weak-coupling limit, and the reduced 
dynamics it generates for the finite-time evolution of the atom. In Sec. Ill, we calculate the 
radiative energy shift of the ground state and the resulting force on the atom. Finally, we 
conclude in Section IV. 



II. THE BASIC FORMALISM 

We consider the evolution in the proper time of a static two-level atom interacting with 
vacuum massless scalar fields outside a Schwarzschild black hole and assume the combined 
system (atom + external vacuum fields) to be initially prepared in a factorized state, with 
the atom held static in the exterior region of the black hole and the fields in the Boulware or 
Unruh vacuum state. When the quantum system of a two-level atom interacts weakly with 
the environment, the reduced dynamics can be obtained by eliminating the environment 



degrees of 
equation 



reedom, yielding an evolution equation of the atom that satisfies the master 



111 ] . We take the total Hamiltonian for the complete system to have the form 

H = H S + H4, + H I , (1) 



where H$ denotes the free Hamiltonian of the atom and that of the environment (a bath 
of fluctuating quantum fields). In fact, can be chosen as the standard Hamiltonian of 
massless, free scalar fields, the details of which are not relevant here. In order to make our 
discussion generic, we will postpone the specification of Hs until later and suppose that the 
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interaction Hamiltonian has the general form 

Ei = p A a <8> B a , (2) 

a 

where A a and B a represent respectively the dynamical variables of the atom and of the 
environment. It should be pointed out that the coupling constant /i should be small, and 
this is required by our assumption that the interaction of the atom with the scalar fields is 
weak. 

Initially, the complete system is described by the total density p to t = p(0) <E> Pb , where 
p(0) is the initial reduced density matrix of the atom, and ps is the state of the environment. 
In the frame of the atom, the evolution in the proper time r of the total density matrix p tot 
of the complete system satisfies 

dr = -iL H [p tot {T)} , (3) 

where the symbol Lh represents the Liouville operator associated with E 

Lh[S] = [E, S] . (4) 

The dynamics of the atom can be obtained by tracing over the field degrees of freedom, that 
is, by applying the trace projection to the total density matrix p(r) = Tr$[p tot (r)] . 

In the limit of weak-coupling which we assume in the present paper, the reduced density 
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is found to obey the master equation in the interaction picture [6 

^ Ps (t) = -i[E LS , p s (r)} + V{p s {r)) , (5) 

where E^s provides a Hamiltonian contribution to the dynamics of the system and is the 
so-called Lamb shift Hamiltonian, whereas, T>(ps(r)) is called the dissipator of the master 
equation. Let us note that here Els is given by 

R ls = P 2 E E SapMAZMM") ( 6 ) 

u) a, 

and the system operator can be decomposed as {f| 



with 



V ' —U=U) 



where the operator IT(z/) denotes the projection onto the eigenspace belonging to the eigen- 
value v of H s . It is easy to see that 

[H a ,At(u>)Ap(u>)] = 0, (9) 
that is, Hls has the same eigenstates as Hs- Meanwhile, the function S a p(u) can be written 



as js] 

S a p{uj) = ~Gap{u) - iT a p(u) , (10) 
where G a p(oo) is the Fourier transform of the reservoir correlation function ((B^(s)Bp(O))) 

/GO 
dse^ s (B+(s)Bp{0)} (11) 
■oo 

and r Q/ g(c<;) denotes the one-side Fourier transform 

roc 

T a p(u) = / dse^ s (B+(s)Bp(0)) . (12) 
Jo 

Then, with the help of 

—J— = P- ± in5(x) , (13) 
1 =F «e x 

where P denotes the Cauchy principal value, it is easy to prove that 



(15) 



The dissipator of the master equation takes the form 

V(p s ) =» 2 J2I2 M") (^(w)ps^(w) - i{A+(a;)^(a;), p s } 
The Hamiltonian H g of a two-level atom can be generically written as [l0| 

3 

2 i=i 
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where Oi (i = 1, 2, 3) are the Pauli matrices, ujq is the energy- level spacing and n = (m, n 3 ) 
is a unit vector. Now, we let n = (0, 0, 1) (i.e., Hs = uoas/2 ) for simplicity. The interaction 
Hamiltonian Hj is taken as 

3 

iZj = /i^<7 a $ a (x) , (17) 

where o"o is the unit matrix and external fields are represented by $> a (x). Comparing it with 
the general form Eq. (J2J), we find that A a = a a , B a = $ Q (x) . 

For the sake of simplicity, we now assume that the field correlation functions are diagonal 
such that 

G + {x-y)6 a p = (* a (x)$p(y)) . (18) 
Therefore, its Fourier and Hilbert transforms can be expressed respectively as 

/oo 
dre^G + (x(r))=5 a ^ a(3 (u), (19) 
-oo 

and 

K.{u)) = — I d\^- = 2i5 a pS af3 {u) . (20) 

m J-oo A - w 

It should be pointed out that here the repeated index does not mean summation. Letting 
|+) and |— ) denote the excited and the ground states, respectively, we have 

A a (u ) = A+(-u ) = \-)(-K\+)(+\ . (21) 

Since the energy spacing is oj , according to Eq. (jSJ), the summation over u in Eq. (jSJ) just 
contains two terms; namely, u = Uq and u> = —ujq. As a result, Eq. ([6]) can be rewritten as 



2 3 



tLjc M j2\+)\{+K\-) (+|-^/c(-c )^|-) <-K|+> 
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.1irt. . M i \ / , I „-..2i 



= -i^K{u Q )\+){+\-i^K{-u Q )\-)(-\ (22) 

This shows that the energy-level shifts of the ground state and excited state are SE^ = 
— ifi 2 JC(— loq) and 5E + = —ifi 2 JC(uj ), respectively, which are caused by the coupling of the 
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atom to the vacuum fluctuation of the fields. The relative energy shift (the Lamb shift) is 
then 

A = ifi 2 [JC(-uj ) - K{uj Q )\ . (23) 



III. CASIMIR-POLDER-LIKE FORCE FOR BOULWARE VACUUM 



In this section, we apply the open quantum system formalism developed in the preceding 
section to address the issue of finite-time evolution of a static two-level atom interacting with 
vacuum scalar fields outside a spherically symmetric black hole and calculate the energy-level 
shift of the ground state. The line element for a Schwarzschild black hole reads 

2M\ dr 2 



ds z 



1 - 



- \dr 



- r\dB l + sin" 



(24) 



r J 1 - 2M/r 

Let us start with the Boulware vacuum which is deemed to be the natural vacuum outside 
a massive body and reduces to the usual Minkowski vacuum at infinity. The correlation 
function for the scalar field in the Boulware vacuum is 



1 1 II 



oo g-iuAt 



Attuj 



Y lm (0A)\ 2 [\i(co,r)\ 2 + \%(co,r)\ 2 ]dLO 



(25) 



and the corresponding Fourier transform with the proper time r of the atom reads 



e tXT G B + [x(r)}dT 



£ 

ml 



du 
2uj 



5(\-u/ v / l-2M/r)\Y lm ( 



; 0^,r)| 2 +|:fe(u;,r)| 2 ] .(26) 



Using Eq. ( 1201) . we can calculate 5E_ for the Boulware vacuum; for convenience, we denote 
it by SE^. Since the exact form of functions Ri(u,r) is not known, a generic expression for 
5E B is not possible to obtain. In what follows, we only examine the behavior of SE B at two 
asymptotic regions, i.e., close to the black hole horizon and at infinity. In order to do this, 
let us recall that [14 1 

Au 2 



oo 

^2(2l + l)\%(u,: 



I _ 2M ' 



1=0 



-. oo 

-;>>/ + 1) i 



r 2M 



r — y oo 



(27) 



1=0 



{ oo 
-V(2/+1)| J B,H| 2 , r^2M, 
r 1=0 (28) 
4w 2 /(l - 2M/r) ~ 4w 2 , r -»■ oo , 

Here Bg is the transmission amplitude. The energy-level shift SE B in two special cases can 
then be calculated. We find, both when r — > 2M and r — > oo, that 



5P/ 8 = -iffJCut-uo) w 5£ + 5^0r • (29) 



Here, we have defined 



5E = -V— —~d\, (30) 

27T^ J A + Uo 



a 2 P F°° A 
5E Qr = -?— f(\,r)-—d\, (31) 
2^ y A + u 



where 



f(X,r) = ^ 1+2,) [ B '^'' , (32) 



4AV 2 



and P denotes the principal value. For the purpose of estimating 5E_ approximately, we 
use the geometrical optics approximation to evaluate the transmission amplitude B^p) jl5f ] . 
In this approximation, it is easy to verify that, if Mp ^> 1, then transmission ceases when I 



exceeds yTJMp and this holds even when Mp is small 15]. So we have 

B t {p) ~ 6(V2JMp - I) , ^(1 + 2l)\B e (p)\ 2 « 27MV • (33) 

This leads to 

Therefore, we can rewrite SE B in both two special cases as 

5E B * 6E + ^^5E . (35) 

Here, 8Eq is just the shift of the ground state energy level in the Minkowski vacuum[2J. So, 
5EE contains two divergent terms, both of which contain a linearly divergent factor SEq. 



According to the Bethe's mass renormalization method 16(, SEq can be written in two parts 



6E = SE' + 6E" = -ff r d\ + ^ r Y^rfA , (36) 

where 8E'q is the energy of a free electron due to its coupling to the field, which can also 
be interpreted as the element stemming from the renormalization of the mass in the kinetic 
energy of the system Hamiltonian. However, it will not contribute to the observed level shift. 
The second part 

^"0 = f tt^A , (37) 

although logarithmically divergent, will give an observable contribution to the atomic level 

the upper limit of integration. For the 



161 ] which was first used in a flat spacetime 



shift after a regularization by taking a cutoff o: 
regularization, we assume that the Bethe method 
can also be applied in the present case. According to Bethe, this cutoff should be taken as 
the electron mass m 1 . Therefore, the renormalized energy- level shift SE B reads 



6E' 



A ^0[>^] (38) 



2ir 2 Vwo- 

It is easy to see that the regular finite energy-level shift is position-dependent, and so gener- 
ates a force on the atom besides the gravitational force. This is a force which has a quantum 
origin and it is in fact a result of the modified vacuum fluctuations due to the spacetime 
curvature. This position-dependent energy shift gives rise to a force on the atom which can 
be calculated by taking the first derivative with respect to r. Therefore, we can obtain the 
Casimir- Polder-like force on an atom outside a spherical star in two asymptotic regions (i.e., 
r — > 00 and r — > 2M): 

^ B d(5E B ) _ 27 / u 2 M 2 u; 0/ „ o*^,„( m \ 



(r - 3M) In — (39) 



dr Ait 2 r A va; / 

A physical realization of the Boulware vacuum would be the vacuum state outside a massive 
spherical body of radius only slightly larger than its Schwarzschild radius. According to 



1 It is worth pointing out that, in flat spacetime, the same result as Bethe's can be obtained if one employs 
a fully relativistic quantum field theoretic approach where no cutoff is present. See Refs. 17 - 2o| 
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Eq. ( )39|) . the Casimir- Polder-like force outside a massive spherical star can be either attrac- 
tive or repulsive. In fact, the force is attractive close to the horizon and repulsive at the 
spatial asymptotic region with a behavior of 1/r 3 . The turning point happens near r « 3M 
where the vacuum field modes are scattered the most. It is interesting to note that r = 3M 
is the location of the unstable circular orbit of photons 



IV. CASIMIR-POLDER-LIKE FORCE FOR UNRUH VACUUM 



For a static atom outside a black hole in interaction with fluctuating massless scalar fields 
in the Unruh vacuum, let us note that the field correlation function is given by [13I. fl4| 



Gu + (x,x') = ^2 



00 g-iuAt 



Attoj 



Y, 



lm\ 



|i^,r)| 2 



6(u)\%{u,r) 



(40) 



where k = 1/(4M) is the surface gravity of the black hole. The corresponding Fourier 
transform with the proper time r of the two-level system reads 



e tXT Gu + {x,x')dr 



1 00 

— y 

1=0 



f(Av^oo)(l + 2Z)| Ri{\y/goo,r) \ +■ 



I _ e -2 7 rA v ^J/K 



• (41) 



A similar calculation as that in the Boulware vacuum yields the radiative energy-level shift 
of the ground-state in two special cases: 



5E U _ = -in 2 )C v (-u) ) 



6E + 6E T + 6E 0r , (r ~ 2M) ; 
SE + SE Tr + SE 0r , (r -> 00) . 



(42) 



where 



and 



5Ei 



fi 2 P f°° A , 1 
: v" 



2tt 2 ./ n 1 - e A / T v A - u X + u 



SE- 



Tr 



tfP [°° /(A,r) , X 



2tt 2 J l-e x / TK X-uj X + loq 



X 



)dX , 



-)dX 



(43) 



(44) 
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Here, T is given by 



K .■>■. 



T) 



H 



•> , ^- ^ ( 45 ) 

with Tf/ = k/2tt, being the usual Hawking temperature of the black hole. This equation is 



the well-known Tolman relation 



21 



22] which gives an effective temperature as measured by 



a local observer. A few comments are in order now. First, 5Et, which is a contribution to 
the energy-level shift of the atom close to the horizon, is structurally similar to the energy- 
level shift in a thermal heat bath at temperature T, while SEtt, which is a contribution to 
that at infinity, shows the effect of backscattering of the thermal radiation off the spacetime 
curvature, represented by the grey-body factor f(X,r). These two terms clearly support the 
notion that there is a thermal radiation flux emanating from the black hole event horizon. 
Second, the other two terms, (i.e., 5E and SE 0r ) are just terms contained in SE^. 

To analyze the force on the atom in more detail, let us note that 5Et can be approximated 
in the low- and high-temperature limits as follows 



5E, 



2vr 2 



A 



1 



A / T v A — u A + u 



tfT 1 /iVT 4 



6cu 
2tt 2 



(too > T) ; 
(w < T) ■ 



(46) 



According the Eq. (131)1 . we can estimate 5Et t in the same limits 

9fi 2 M 2 T 2 g 00 9/i 2 vr 2 M 2 T%o 



5Exr = — ^°° SEt 
4r 2 



3Co>0 r 
87r 2 r 2 11 \ T 2 



20u$r 2 



("o > T) ; 



(w < T) . 



(47) 



Since close to the horizon, g 00 approaches zero, T = Tn/y/gm 3> wo is always satisfied. 
Consequently, the position-dependent energy-level shift can be written as 



5E 



u 



2tt 2 



In 



u; 

9007^2 
H 



fi 2 u i in 
+ „ m 



2tt 2 



w 



1 + 



27M 2 (7oq 
4r 2 



(4£ 



When the atom is held static in the spatial asymptotic region, i.e., when r — > oo, T ~ Tjj . 
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Now, the energy-level shift becomes 



5E U 



9fi 2 M 2 T 2 l_ + f^o ln (m^ 



8w 



27/i 2 M 2 w 
87r 2 r 2 



2tt 2 



1 + 



27M 2 g 00 
4r 2 



#00 hi 



CO, 



2 -, 



9oo7pr 
1 H 



/rw . (m 



j ) 



1 + 



27M 2 g 00 
4r 2 



(w > T H ) ; 



(w « Tff) . 



(49) 



This position-dependent energy shift gives rise to a force on the atom which can be calculated 
by taking the first derivative with respect to r. Close to the horizon, we find that 

d(6E") fi 2 Mu 27/i 2 w 



F u = 



1 (™\ 
64Mtt 2 n KuqJ ' 



(50) 



9r n 2 (r-2M)r 
So the Casimir-Polder-like force is attractive and actually diverges at the event horizon. Let 
us note that the classical force that is needed to hold the atom static at the horizon also 
diverges. If the atom is in the spatial asymptotic region, this force can be approximated as 



d(SEE) 
dr 



9fi 2 M 2 T H 2 1 27fi 2 M 2 u 1 



4Wn 



+ 



4tt 2 



In 



/ m 



\u 



)■ 



27 fi 2 M 2 u 1 (T H 



2ix 2 



In 



u 



| 27 fi 2 M 2 u 1 h , 



4tt 2 



/m \ 

\UJ J 



(w < T H ) . 

(51) 



Typically, one has m ^> Wo, so, far from the black hole, the force will be attractive if Th ^> m, 
and repulsive otherwise. It is interesting to note that the contribution to the force due to 
the presence of the Hawking radiation is always attractive. Therefore, the collapsing of a 
massive star into a black hole and the thermal radiation generated as a result makes the 
Casimir-Polder-like force more attractive than repulsive. 



V. CONCLUSION 



In summary, we have calculated the energy-level shift of a two-level atom outside a spher- 
ically symmetric black hole in the paradigm of open quantum systems by looking at the 
time evolution of the atom interacting with massless scalar fields in the Boulware vacuum 
and the Unruh vacuum. The time evolution of the atom is governed by a master equation 

12 



obtained by tracing over the field degrees of freedom from the complete system. Our results 
show that, for an atom in the ground state, the level shift is position-dependent and gives 
rise to a force on the atom besides the classical gravitational force. 

For the case of the Boulware vacuum, which represents a star that has not collapsed 
through its event horizon, this force is attractive near the horizon, and is repulsive far away 
from the black hole with a behavior of r -3 . The turning point occurs near r ~ 3M where 
the vacuum field modes are scattered the most. For the case of the Unruh vacuum which 
represents a radiating black hole, we find that the contribution to the Casimir-Polder-like 
force due to the presence of Hawking radiation is always attractive, and in fact this attractive 
force diverges as the horizon is approached. 
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